TeV Scale Singlet Dark Matter 



Eduardo Ponton" and Lisa Randall 6 , 

a Department of Physics, Columbia University, 
538 W. 120th St, New York, NY 10027, USA 
^Jefferson Laboratory of Physics, Harvard University, 
Cambridge, Massachusetts 02138, USA 



Abstract 

It is well known that stable weak scale particles are viable dark matter candidates 
since the annihilation cross section is naturally about the right magnitude to leave the 
correct thermal residual abundance. Many dark matter searches have focused on relatively 
light dark matter consistent with weak couplings to the Standard Model. However, in 
a strongly coupled theory, or even if the coupling is just a few times bigger than the 
Standard Model couplings, dark matter can have TeV-scale mass with the correct thermal 
relic abundance. Here we consider neutral TeV- mass scalar dark matter, its necessary 
interactions, and potential signals. We consider signals both with and without higher- 
dimension operators generated by strong coupling at the TeV scale, as might happen for 
example in an RS scenario. We find some potential for detection in high energy photons 
that depends on the dark matter distribution. Detection in positrons at lower energies, 
such as those PAMELA probes, would be difficult though a higher energy positron signal 
could in principle be detectable over background. However, a light dark matter particle 
with higher-dimensional interactions consistent with a TeV cutoff can in principle match 
PAMELA data. 







1 Introduction 



Dark matter has received a lot of attention of late as new dark matter searches ramp up. Of 
particular interest is the increasing capacity to detect dark matter in both direct and indirect 
channels. The latter rely solely on dark matter annihilation, which is nice in that it doesnt 
assume any particular type of interaction with the Standard Model and furthermore the anni- 
hilation rate is generally connected to the annihilation cross section responsible for the current 
dark matter abundance. 

Given the importance of dark matter searches and our lack of knowledge as to the true 
nature of dark matter, it makes sense to explore the range of possibilities and what their 
implications would be for current and future detectors. In this paper we will consider singlet 
dark matter candidates with mass of order one to a few TeV. We assume a Z2 symmetry 
that prevents any operator allowing decay and therefore ensuring stability. This is perhaps the 
simplest dark matter candidate there can be. In fact, such a possibility has been previously 
considered in Ref. [H El El Hj , but in a lower mass region. In this paper we concentrate on the 
remaining allowed mass range, of order one to a few TeV, which phenomenologically is also a 
viable possibility. We concentrate on some novel scenarios that arise in a framework with a low 
cutoff scale. 

Although we mostly take an agnostic approach about the source of this dark matter, we 
also focus on TeV scale particles that arise in a theory with a TeV cut-off scale. Such a scenario 
can occur for example in the RS framework [5J. See also [6] for a study of heavy DM in a 
supersymmetric theory with a relatively low cutoff. 

In this paper we show the range of allowed parameters giving the right relic density and then 
consider whether such dark matter has any chance of being detected. We find that annihilation 
into photons might provide a visible signal at high energy gamma ray detectors such as HESS 
or VERITAS, particularly if higher- dimension operators are present. We also consider more 
model-dependent scenarios in which annihilation into positrons can also occur. We show the 
signal can exceed background with reasonable assumptions, but most likely not in the PAMELA 
range for a TeV-scale mass. 

On the other hand, we observe that a dark matter candidate of about 100 GeV whose dom- 
inant decay mode involves direct positron emission, such as can occur with a higher- dimension 
operator suppressed by the TeV scale, matches PAMELA data quite nicely. 
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2 Singlet Dark Matter 



We start by discussing the relic density computation for a thermally produced Standard Model 
singlet. We consider first a renormalizable four- dimensional theory This analysis would of 
course also apply to a nonrenormalizable theory so long as the renormalizable coupling of the 
singlet scalar to a Higgs dominates annihilation, including a five- dimensional theory with a 
brane-bound scalar or any five-dimensional theory where the higher- dimension operators are 
suppressed. 

We then consider a more exotic possibility that could in principle give rise to a detectable 
positron signal. We will see this scenario is unlikely to explain the PAMELA data, although it 
could give rise to a detectable signal in the high-energy positron range. 

2.1 Thermal Relic Abundance for a Singlet 

We assume a singlet field $ protected by a discrete Z 2 symmetry $ — > — $ in a nonrenormaliz- 
able theory with a TeV cutoff scale, A. Without any additional fields, the only renormalizable 
operator that involves SM fields is 



where H is the Higgs doublet and A is a dimensionless coupling. Such an operator can arise in 
an RS scenario for either IR brane- localized or bulk scalars $. For an IR brane-localized scalar 
(assuming the Higgs is also IR localized), the corresponding operator is 



For a bulk scalar <3> the operator Eq. (CEJ) can be induced from a non-renormalizable operator (to 
be discussed in the following subsection). If the cutoff is at the TeV scale the effective coupling 
A can easily be of order one, so that the following analysis applies. 

The interaction in Eq. ([1]) can lead to the direct self-annihilation of $ particles into a pair 
of Higgses, and also, if the annihilations occur after the electroweak phase transition, into pairs 
of SM gauge bosons and fermions through s-channel Higgs exchange. When the $ mass is 
much larger than the Higgs mass, the direct annihilation into Higgses dominates (this includes 
annihilation into the Goldstone modes, hence the WlWl and ZlZl channels). Annihilation into 
two Higgses in the limit that M$ ^ v ew (with i>ew the Higgs VEV) gives in the non-relativistic 




(1) 




(2) 
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regime 

(VI^HHV) - ^J, (3) 

where v is the relative velocity of the annihilating particles, and the brackets denote thermal 
averaging. 

Notice that other annihilation channels have to proceed through operators suppressed either 
by the cutoff scale A or by a loop factor. As we will argue in the next subsection, when the 
$ particles propagate in the bulk of an RS scenario those channels might be relevant (and 
could even dominate depending on couplings) in the total self-annihilation cross section, and 
therefore in the determination of the relic density. However, for a conventional four-dimensional 
scalar (or for a brane- localized $ in an RS scenaric*]) all other channels are expected to give 
a relatively small contribution when the $ mass is less than the cutoff scale. For instance, 
annihilation into SM fermions would proceed through operators that also involve the Higgs 
field, of the form Q 2 Hijj 1 'ip2, and are suppressed at least by order (vpw/^) 2 , where vew is the 
Higgs vacuum expectation value and A = Ae~ kL is the warped down cutoff scale if in a 5D 
warped framework, or more generally the cutoff scale of the 4D theory. Decays into SM gauge 
bosons are also expected to be subdominant so long as the DM mass is less than the cutoff, 
even if the cutoff scale is low, and will be discussed in Subsection 13.11 in the context of DM 
indirect signals [see Eqs. ([TBI and (TT4l) and ensuing discussion]. 

Under the assumption that the DM candidate is heavy (say 1 TeV or so) and is thermally 
produced, the DM relic abundance is controlled by Eq. ([3]). Taking into account only the 
annihilation into Higgses through the operator Eq. (JJJ), and requiring that the observed DM 
abundance is completely accounted for by $ particles, we can determine the coupling A as a 
function of M$ from the WMAP constraint Qp>Mh 2 ~ 0.11 [7j and 

„ ,„ 1.04 x 10 9 GeV" 1 x F 1 
M P Jgl {av) 



where Mp rs 1.22 x 10 19 GeV is the Planck mass, xp = M^/Tp, with Tp the freeze-out 
temperature, g* is the effective number of relativistic degrees of freedom at freeze-out, and 
(av) is the thermally averaged annihilation cross section times relative velocity in units of 
GeV -2 . We have assumed that the $ particles are thermally produced and remain in thermal 



-"■The operators discussed in Subsection 12.21 vanish for a brane localized <£> since always one of the chiralities 
of any bulk fcrmion '5 vanishes on the brane. 
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Figure 1: Left panel: annihilation cross section, (a^^HH'v) , for a brane-localized scalar in the 
non-relativistic regime as a function of M$, imposing the WMAP constraint on the DM relic 
density. The arrows indicate the points where the freeze-out temperature (~ M$/25) crosses 
the W ± and Z° thresholds. Right panel: the corresponding coupling A, defined in Eq. ([1]), as 
a function of M$. 

equilibrium until freeze-out, which requires a coupling A > 10 -8 [3]. For masses M$ in the 

few TeV range, A is always of order unity, as shown in the right panel of Fig. [TJ so the above 

assumption is self-consistently satisfied. 

The required (non-relativistic) cross section as determined from the WMAP constraint to 

be (av/c) ~ 0.8 pb and is shown in the left panel of Fig. [1] as a function of M$, where the 

(weak) dependence of x_p on the cross section and the effective number of relativistic degrees of 

freedom is included^ Throughout the range of interest we have xf ~ 25, while g* is of order 

90. The above simple picture is rather generic for a stable TeV scale scalar field whenever the 

effects of higher-dimension operators can be neglected. However, when the scale suppressing 

the higher-dimension operators is near the TeV scale, other more exotic scenarios are possible. 

Such a situation, though less likely, could arise within the RS framework, and will be illustrated 

with a bulk scalar in the next subsection. 

2 The conversion factor from GeV -2 to pb is 0.3894 x 10 9 GeV 2 pb. Also, to convert the above cross section 
from pb into units of cm 3 s~ 1 [the CGS units for (av)] one should multiply the number in pb by (10 -36 cm 2 ) c ~ 
3 x 1CT 26 cm 3 s- 1 . 
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2.2 Bulk Singlet Dark Matter 



Another natural possibility in a 5D warped background is that the DM candidate arises as the 
lightest KK mode of a bulk scalar. In order to be concrete, and simply for illustration purposes, 
we will assume in the following that there is a bulk SM singlet scalar obeying (— , +) boundary 
conditions (Dirichlet on the UV brane, Neumann on the IR brane). In this case, the mass of 
the lightest KK mode is determined by only two dimensionless parameters, and can be easily 
below those of the gauge KK resonances (say around 1 TeV), as discussed in more detail in 
Appendix |A] We also assume that the SM fermions and gauge fields arise from bulk fields. 

The couplings to (an IR localized) Higgs field proceed now through the higher- dimension 
operator 



where A is the cutoff scale and A' is a dimensionless coupling. After KK reduction, this induces 
a coupling of $s to Higgses similar to the one discussed in the previous subsection with the 
identification A = A'/| (1) /(AL) « X'(2k/A). Here / $ (i) « \JlkL is the t&W wavefunction 
evaluated on the IR brane, where is the lightest scalar KK mode (the DM candidate). If this 
channel dominates the self-annihilation cross section, the relic density computation proceeds 
in exactly the same way as in the case of a brane-localized scalar discussed in the previous 
subsection. As was mentioned there, the observed relic abundance requires an effective 4D 
coupling A of order one. Notice that for a bulk scalar, in spite of the suppression k/A, this is 
easily consistent with the NDA bound A' < 24tt 3 [8]; in fact, for k/A ~ 1/10 the fundamental 
coupling A' is well into the perturbative regime, so that the computation is under theoretical 
control. 

It is possible, however, that channels other than the annihilation into Higgses are important 
or even dominate, which could in principle differentiate between brane and bulk dark matter 
candidates. This scenario requires the value of A' well below its NDA value with other couplings 
closer to what NDA would suggest. Since we are taking an agnostic attitude and are interested 
primarily in potential signatures and ways to identify the various possible scenarios, we consider 
this possibility next. 

For this analysis it is useful to rewrite the annihilation cross section into Higgses as 



(5) 




5 



where A = A e~ kL is the warped down cutoff (of order a TeV) and similarly for k = k e~ kL . Now 
consider operators involving a 5D fermion field (giving rise to a SM fermion as its zero-mode), 
for instance 



$ ^>\I> (7) 

2A 2 ' 1 ' 

where \I/ is the bulk fermion and is a dimensionless coupling. This operator leads to the 
annihilation of $ particles into a SM fermion and one of its KK resonances, e.g. 

^-($«) 2 fe (1 V (0) +? (0 V (1) 

2A(AL) V I ■ 

where ip^ is the first KK mode of the bulk fermion and ip^ is its zero mode (with a well- 
defined chirality). The effective 4D coupling depends on the various extra-dimensional profiles 
through 

V = ^J L dye k ^fl w f^ } f^ , (9) 

where all the wavefunctions are normalized as in Eq. (1321) of Appendix lAl 

Assuming that the channel — > tp^ip is open, i.e. Mmi) +m^ ) < 2M$, the correspond- 
ing annihilation cross section is |f 

a ^ { ° )V = 167rA 2 (AL) 2 M^l 2 ' 

where N c = 3 for quarks while N c — 1 for leptons and, for simplicity, we neglected the zero-mode 
mass 77ty(o). In the non-relativistic limit one has a^w^v = a + bv 2 + ■ ■ ■ , with 

N c Xl V 2 _ N c \%rf 

8vrA 2 (AL) 2 1 V> ' 64ttA 2 (AL) 2 1 V V ' ' 1 J 

where y = M 2 m /4M$. The related process a^^mv is given by the same expression. 

The magnitude of the annihilation cross sections into fermions depends strongly on the 
localization of the fermion zero- mode through the parameter r\ of Eq. (Q. Recall that the 
fermion zero-mode wavefunctions are proportional to e ( - 1 / 2 ~ Cf * )ky , where c/ parametrizes the 5D 
fermion mass in units of the AdS curvature scale k. The massive KK mode wavefunctions are 
all strongly localized near the IR brane. There are a number of distinct scenarios according to 
how flavor is generated: 



3 To simplify notation we will refer to the DM candidate simply as <£>. 
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1. The SM fermion mass hierarchies arise from the exponential wavefunction localization and 
the overlap with an IR localized Higgs field. This scenario has the advantage that both 
calculable flavor changing effects (from KK gluon exchange), as well as non-calculable 
effects from flavor changing non-renormalizable operators, are significantly suppressed [91 
[TUl [TTl [12] . One expects the third generation quarks (most likely the right-handed top) 
to couple most strongly to $. 

2. All fermions share the same parameter Cf, and are localized close to the IR brane 
(cf < 1/2), so that their couplings to $ are sizable. Somewhat more generally, EW 
precision constraints allow different localization parameters for different fermions so long 
as those fermions having identical quantum numbers have nearly the same c/ (when IR 
localized; otherwise we are in scenario 1 above). In these scenarios, as-yet unspecified 
flavor-violating interactions would be necessary to explain the fermion mass hierarchies, 
while not generating dangerous FCNC effects from higher- dimension operators suppressed 
by the TeV scale. 

3. Fermion mass hierarchies arise from localization in the extra dimension but the Higgs 
field is located on or near the UV brane (for example, if the Higgs mass is stabilized 
by supersymmetry (SUSY) and SUSY breaking is connected to the IR scale). In this 
case, the lightest fermions would be localized closer to the IR brane and have the largest 
couplings to $. 

Among the fermion channels, $ annihilates dominantly into the fermions closest to the IR 
brane, since the $ wavefunction is localized near the IR brane. To calculate the annihilation 
rate, we need to estimate the expected size of these couplings, which can then be compared 
to the couplings to Higgses discussed above, or to the annihilation into gauge bosons (see 
Subsection 13.11) . 

For a fermion localized near the IR brane (localization parameter Cf < 1/2, but not very close 
to 1/2), one finds rj ~ (l/5kL)(\ / 2kL) 3 ^(1 — 1cf)kL « \J\~ c f wnere eacri KK wave- 
function contributes a factor \/2kL, the last factor corresponds to the c/-dependent zero-mode 
wavefunction, and the factor l/(5kL) is a measure of the region that contributes to the integral 
in Eq. ©3 Compared to the annihilation into a pair of IR localized Higgses, Eq. ([6]), the annihi- 
lation into fermion and KK fermion is "suppressed" by order (N c /2)(\^/\') 2 (M$>/A) 2 (l/AL) 2 , 

The factor of 1/5 is determined by comparison to the exact result, Eq. ([5]), and reproduces it within 30% 
for -0.5 < c f < 0.4. 
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Figure 2: Annihilation cross section, (a^^m^mv), for bulk DM, as a function of M$, imposing 
the WMAP constraint on the DM relic density. The curves marked as "freeze-out" correspond 
to the annihilation cross section at the time of freeze-out (where the typical velocities were 
of order v/c ~ a/2/25 ~ 0.3). The lower curves correspond to the annihilation cross section 
in the ultra non-relativistic regime, as would be relevant for today's conditions. The various 
curves correspond to different choices of the fermion localization parameter Cf that controls 
their masses and couplings. The arrows indicate the points where the freeze-out temperature 
(~ M$/25) crosses the W ± and Z° thresholds. The curves are terminated (with black dots) 
when = 247T 3 , which we define as the strong coupling regime (see text). We assume that 
A = 8k. 

where it was assumed that the annihilation into fermions is not near threshold, and we used our 
estimate for rj and take Cf of order one. The NDA estimate for is 24-7T 3 , which is the same 
as for A'. However, the discussion after Eq. (jSJ) indicates that a correct thermal relic abun- 
dance requires a much smaller coupling A' < A/ (2k). Taking A < 10k, M$ ~ k, kL »s 34, 
one can see that the annihilation into an IR localized fermion and its lightest KK mode could 
dominate the annihilation cross section of $ particles. Other operators that contribute to the 
self-annihilation cross section are expected to give subdominant contributions when M$ <C A. 
Nevertheless, the operators that lead to annihilation into gauge bosons can be interesting from 
the point of view of DM signals, and are discussed in subsequent sections. 

We explore here the exotic picture where the $s are kept in thermal equilibrium dominantly 
by annihilation into a fermion and its lightest KK mode. For concreteness, we imagine here 
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scenario 1 discussed above, where the annihilation proceeds mainly into tops and its lightest 
KK resonance, but the same would hold in scenario 3 with one of the lightest leptons (either 
electrons or neutrinos) replacing the top. The results can also be applied in a straightforward 
way to scenario 2 with all fermions localized identically: one should just include a multiplicity 
factor 3 x [3 x 4 + 3] = 45. 

In scenario 1, the processes taken into account are 

$$ _^ T^k R and -> T { L 1] t R , where 
T"W is the first KK excitation of the RH top tower. The annihilation cross section depends on 
the lightest KK scalar and fermion masses M$ and M T (i>, which are both of order k. We fix k 
and obtain different values of M$ as described in Appendix [A] The KK fermion mass has some 
dependence on q, which controls the localization of the t R wavefunction. The overall strength 
of the cross section depends on the combination X t /[A(AL)] = \ t /[k{kL)}{k/ 'A) 2 . Assuming 
again that the $s account completely for the observed DM energy density, one can then fix 
the quantitj^l \ t {k/K) 2 using the WMAP result and Eq. (TjJ with (av) = a + 36/ Xf where 
the coefficients a and b are given in Eq. (flTT) . In Fig. [2] we show the result for several values 
of the fermion localization parameter Cf = c t (for the RH top most likely c t is close to 0). 
As expected, the annihilation cross section at freeze-out is (av/c) ~ 0.8 pb. However, unlike 
the case of annihilation into scalar particles such as the Higgs field discussed in the previous 
subsection, both the a and b terms give a comparable contribution. As a result, the annihilation 
cross section at very low-temperatures, being dominated by the a term, is somewhat different 
from the cross section at freeze-out. This is relevant for annihilation under today's conditions, 
and is also shown in Fig. [2j The curves marked as "At freeze-out" correspond to the annihilation 
cross section at the time of $ decoupling (when the typical velocities were v/c ~ a/2/25 ~ 0.3), 
while the curves in the lower part of the plot correspond to the annihilation cross section in the 
ultra non-relativistic regime, and correspond essentially to the a-term in Eq. ( fill) . The behavior 
observed in these curves arises from the fact that as AM = 2M$ — M tW approaches zero, the 
annihilation cross section vanishes. Specifically a ~ A 2 (AM) 2 and b ~ A 2 AM. Thus, near 
threshold the b term dominates, and the WMAP constraint requires the scaling A 2 ~ 1/AM. 
This explains why the annihilation cross section at very low temperatures decreases as M$ 
decreases (for fixed M Tm ), since a ~ A 2 (AM) 2 ~ 1/A 2 ^ ~ AM. 

We terminate the curves at the point where the coupling A^, reaches the strong coupling 
value given by NDA, A^, ~ 24-7T 3 , assuming A ~ 8k. As explained above this happens near the 

5 If A is defined as the scale where the SU(3)c gauge factor gets strong, then NDA gives AL ~ 247r 3 / (3g^) , 
where g s is the 4D color coupling at the KK scale. For kL w 34 this corresponds to A/k w 8. 
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threshold for top-KK top production. The different curves are terminated at different points 
due to the c t dependence of the KK fermion mass M Tm . Thus, at strong coupling, A^ cancels 
the volume suppression factor AL in Eq. (flOl) that arises from the fact that the operator Eq. (JTj) 
is suppressed by two powers of A. The fact that this channel then dominates over the Higgs 
pair production channel, in spite of arising from an operator of higher dimensionality can then 
be understood as due to the strong localization near the IR brane of the RH top quark, as 
encoded in the parameter r] of Eq. (Q as well as the different values of the couplings of the 
associated operators. Away from threshold the coupling is a factor of 5-10 below the NDA 
value, so that the perturbative computation can be trusted. 

It is therefore plausible that the annihilation into Higgses plays a sub dominant role in the 
determination of the DM relic density. Of course it is straightforward to take both channels 
into account when they give a comparable contribution, but we will not do so here and turn 
instead to the possible DM signals of these scenarios. Note however that sizable brane-localized 
kinetic terms (that were not included in the above analysis) are known to lower the lightest 
KK masses significantly [13]. Thus, even for k = 1.2 TeV (as is suggested by the EW precision 
constraints as a lower bound on k, and as assumed in Fig. [2]) the KK masses can easily be 
somewhat below a TeV. The qualitative behavior of the ultra non-relativistic cross section 
persists: it is of order 0.8 pb, and decreases by a factor of about two near the threshold for 
$$ — > T^H annihilation (assuming that this is the main annihilation channel and that we are 
in the perturbative regime). Thus, in the following phenomenological analysis, we will allow a 
large range of KK masses and analyze the consequences for indirect detection. 

3 Indirect Detection 

Because the scalar couples to the Higgs, interactions relevant for direct-detection experiments 
are in principle possible [3]. However for the heavy scalars we are talking about the direct 
detection rate will be too low so we concentrate on indirect signals. 

In this section we consider such possible signals from the DM candidates discussed in Sec- 
tion [2j We will argue that current experiments may be sufficiently sensitive to detect photons 
(or possibly positrons) from dark matter annihilation, most likely when non-renormalizable 
operators are present. We will present our bounds in terms of constraints on the cutoff scale A 
appearing in these operators, which in the RS context can be understood as being related to 
the fundamental gravity scale and more generally represents a scale of strong interactions. 
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We first concentrate on the most distinctive signals, — > 7X and — > e + X, where the 
photon(s) and positron are produced from direct 2-body decays and have well-defined energies. 
We also consider the more exotic decay chain involving a KK lepton, which generally yields a 
continuous spectrum (even before propagation thorugh the interstellar medium) except when 
this KK lepton is sufficiently heavy to be produced almost at rest so that the positrons that 
result from its decay have a spectral distribution similar to those of primary positrons. 

Subsequently we will consider possibilities from the decay of the Higgs that would occur as 
a consequence of the dark matter-Higgs coupling. 

3.1 Photons 

We now consider possible photon signals arising from annihilating dark matter. For a continuous 
photon signal the total flux is obtained by integrating from some detector-dependent threshold 
energy up to the DM mass. In the scenarios discussed in Section [2] the continuous signal is 
likely too small to see but we comment on such decays at the end of this subsection. 

We start by discussing the more interesting signal arising from the direct decays of the 
(slowly moving) DM particle into photons proceeding from higher dimension operators, in which 
case the final photon is nearly monoenergetic. Both decays into two photons and a photon and 
a Z could in principle contribute. The photon energy in the first process is approximately 
equal to M$, while in the second process it is approximately M$(l — M§/4M|). For DM in 
the TeV range, the energy resolution of ACT's is not enough to resolve the two lines and they 
both appear to have energy essentially equal to M$. It is therefore appropriate to add the two 
photon signals in the flux. 

The annihilation of a SM singlet into photons can proceed via higher dimension operators 
which we write as 

2 2 / 

_ LI $2 F F r e JL^_ $2 Z , (12) 

8A 2 As w c w V V ' 

where and Z^ u are the photon and Z gauge boson field strengths, sw is the sine of the 
weak mixing angle, A is the effective 4D cutoff scale, and couplings of order one. 

The operators in Eq. f|T2|) are the 4D effective operators induced by bulk or brane-localized 
operators, depending on whether $ arises from a bulk field or is localized on the IR brane. In 
the RS framework the cutoff might be expected to be around the TeV scale and not far from the 
mass of so that the resulting annihilation into photons need not be extremely suppressed. 
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In the ultra non-relativistic limit (DM particle velocities in the galaxy are of order v 
10 _3 c), the interaction terms in Eq. ( fl2l) give rise to the cross sections 



(a 2j v/c) 



M*\ 4 3vraV 



A J M? 



\ 2 / M - \ 4 



1 TeV\ 2 (Mi 



°- 2pb i^r x |K - (13) 



where a is the fine structure constant, and an annihilation cross section into 7Z 

U\ 4 67raV 2 / M 2 Z 



,'M*V 67raV 2 / M 2 
KWc) « ^ 1 - 



A ; ^c^M 2 V 4M| 
o r ^ /l TeV\ 2 /M \ 4 , 2 

Notice that besides the enhancement in the jZ channel due to the gauge coupling (the factor 
l/s^Cw m 5.6), there is an additional factor of 2 difference due to the identical particle nature 
of the final state photons in the 2j channel. This factor is compensated by the explicit factor 
of 2 in Eq. (TlBl) that accounts for the two photons in the final state. 

The rates given in Eqs. (fl3|) and (fl4"I) are small when the $ mass is low due to the strong 
(M$/A) 4 dependence. However, if A is not much above M$, ground based Cherenkov detectors 
can be sensitive to this signal. We will now interpret current bounds in terms of the implications 
for the cutoff scale A. 

The differential photon flux from a direction that forms an angle ip with the galactic plane 

is 

L_ = \ l a v ) — 1 / dip (r) , (15) 

dttdE ^ X 1 dE 47rM| 7 

where r 2 = I 2 + r 2 , — 2lrocosip, with ro ~ 8.5 kpc the distance from the Earth to the galactic 
center. The integration is along the line of sight, dl, and encodes the information about the 
DM distribution, assuming a spherical DM halo of energy density p(r). The particle physics 
input enters through the thermally averaged cross section times relative velocity (for channels 
labeled by i) and the differential photon yield in channel i, dN^/dE where we add the 77 and 
'jZ signals. We have 

2 

J(Att)Att , (16) 



$ 7 = 5.66 x 1(T 12 cnrV 1 



lpb / V 1 P b 



1 TeV 
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Figure 3: Left panel: J(AQ) as a function of AQ for three different halo profiles (taken from 
Ref. [17]). Right panel: J(AQ) x AO as a function of AO. 

where the factor of 2 corresponds to the two photons per decay in the 2j annihilation channel, 
J (AQ) = (1/AQ) J An J(ip)dQ integrates over the angular acceptance of the detector AO, and 
J (if)) is conventionally defined as 

JN,) = \ ( _L_ -] I" dlp 2 (r) . (17) 

K ^ J 8.5 kpc V0.3 GeV/cm 3 J J K ' K J 

The quantity J(AO) depends on the DM halo profile and can vary over several orders of 
magnitude depending on the halo model, when looking towards the galactic center. It has been 
computed for several DM halo models, and as a function of AO in [TJ] . In the left panel of Fig. [3] 
we reproduce J(AQ) as a function of AQ for three different halo profiles: the Moore et. al. 
profile [15] (a rather cuspy profile), the widely used Navarro- Frenk- White (NFW) profile [16j, 
and a smooth isothermal profile [T7] . For reference we also show in the right panel the product 
J (Aft) x AQ as a function of AQ. The angular acceptance AQ depends on the experimental 
setup. 

Several experiments exist that can search for photons from dark matter annihilation. Among 
the particle community, FERMI (formerly known as GLAST) has recently received a great 
deal of attention. FERMI is a satellite-based detector with excellent angular coverage (greater 
than about 2 sr) and fairly good energy resolution (< 10%) [18]. FERMI should have a flux 
sensitivity of order 10~ 10 photons cm -2 s _1 for photon energies between about 20 and 300 GeV 
(with decreasing sensitivity at lower energies and no sensitivity at higher energies). 
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Ground-based detectors, on the other hand, have much smaller angular coverage and worse 
energy resolution. However their flux sensitivity is similar to that of FERMI at around 50 
GeV [H] and rapidly overtakes FERMI's sensitivity, depending on angular coverage, at higher 
energies. From the right panel in Fig. [31 and taking the NFW halo profile, we see that when 
An = 1(T 5 sr, a typical value used in HESS, one has J(AO) xA!)k KT 1 , while for FERMI 
with AVI = 2 one has J(Aft) xAfift! 20. Therefore, if HESS achieves only An = 10~ 5 , it can be 
more sensitive than FERMI to the photon signal we discuss at energies of about 250 GeV [19] 
(near the end of FERMI's sensitivity range). Ground-based detection will be relatively more 
sensitive with a Moore profile and less so with an isothermal profile. 

HESS and VERITAS could reach a larger angular coverage, since their fields of view (5° for 
HESS and 3.5° for VERITAS) correspond to An ~ 10~ 2 sr. If such angular acceptances are 
reached, and assuming an NFW profile so that J(Af2) x An « 3 (see Fig. [3]), these ground- 
based detectors could overtake FERMI's sensitivity even at photon energies of about 100 GeV. 
Clearly the flux sensitivity is better for either of the ground based experiments^] for reasonable 
dark matter masses above about 100 GeV and the determining factor of which is better is likely 
to be the angular resolution^ 

Of course, without knowing the dark matter profile, it makes sense to search in both satellite 
and ground-based experiments at low energies. However, it should be borne in mind that 
most dark matter models predict a monochromatic photon signal only at one loop so indirect 
detection is unlikely to be sufficiently sensitive to this type of signature of standard thermal 
dark matter. Supersymmetric dark matter annihilation into photons offers perhaps the best 
possible loop-suppressed scenario because the loop can be enhanced [211 122] due to a reasonably 
large numerical factor and because for a higgsino dominated neutralino an enhancement in the 
loop diagram due to near degeneracy with an intermediate state can lead to a cross section that 
saturates with l/m^y dependence (rather than suppression by the potentially bigger dark matter 
mass). This signal is potentially observable, however, only for light dark matter candidates for 
which the flux is big (and where FERMI is sensitive). Otherwise the cross section is too small. 

We note that the direct signal we discuss is at higher energies since we assume a heavy 
dark matter candidate and therefore concentrate on ground-based experiments since they have 

6 VERITAS does not always point toward the galactic center, however, so the flux sensitivy in that regime is 
not guranteed 

7 The sensitivity of ACTs to the photon signal at large Af2 could be limited by the cosmic ray background, 
since subtracting the signal from a nearby region can have a significant effect for shallow profiles [20]. However, 
for more peaked profiles such as NFW, this is expected to be at most an order one effect. 
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An 


J (AQ) x AO, 


A [TeV] at M* = 1 TeV 


2(a 2 - r v) + (cr-yzv) [cm 3 s -1 ] 


10" 5 


KT 1 (NFW) 


2 


5.4 x KT 27 (1.8 x 1IT 1 pb) 


1(T 3 


1 (NFW) 


3.5 


5.8 x KT 28 (1.9 x KT 2 pb) 


any 


10 2 (Moore) 


11.3 


5.3 x KT 30 (1.8 x KT 4 pb) 



Table 1: Sensitivity of HESS or VERITAS to the cutoff scale A for representative Afi's (NFW 
and Moore et. al. halo profiles). A DM candidate with mass M$ = 1 TeV annihilating into 
monoenergetic 1 TeV photons is assumed. We assume k — k' = 1 (see text). The last column 
gives the thermally averaged annihilation cross section into photons for the corresponding A 
(and for M$ = 1 TeV). 



better sensitivity. Although the number density of heavy dark matter particles is lower than 
that for lighter dark matter candidates, our prediction is a tree- level effect, albeit through a 
higher-dimension operator, and the cross section can be bigger than typical supersymmetric 
annihilation cross sections [21], which saturate at about 10~ 28 cm~ 3 s~ 4 . For example, in the 
first row of Table [T] we see that for A = 2 TeV and M$ = 1 TeV the annihilation cross section 
is 5.4 x 1CT 27 cm 3 s _1 . 

For instance, for 1 TeV photons, HESS has a flux sensitivity of about 10~ 13 cm~ 2 s~ 4 . Using 
An = 10~ 5 sr and taking the NFW halo profile, we see from the right panel in Fig. [3] that 
J(Afi) xAO^ KT 1 . The expected flux is then $ 7 = 1.6 x KT 12 cm' 2 s^fl TeV/ A) 4 , which 
could be translated into a bound A > 2 TeV (we also assumed k = k' = l)Jj Under the same 
assumptions, for 2.3 TeV photons HESS would put a bound A > M$ ~ 2.3 TeV. On the other 
hand, the sensitivity could in principle be bigger or smaller according to the dark matter profile. 
For example, for the rather peaked Moore et. al. profile, one has J(AO) x An « 10 2 and the 
non-observation of a line at 1 TeV by HESS would correspond to a bound A > 11.3 TeV. This 
is the expected scale for A in several well-motivated scenarios that take into account the EW 
constraints [23l [2il [251 [26] . 

For HESS or VERITAS operating at An = 1CU 3 sr and using again the NFW halo model 
with J (Ail) x An w 1, the expected flux would be $ 7 = 1.6 x KT 11 cm" 2 s _1 (l TeV/M$) 2 x 



8 For M$ = 1 TeV and A = 2 TeV, the non-relativistic annihilation cross section into jZ is (a 1 zv/c) w 
0.15 pb, which is smaller than the cross section necessary for M$ to account for the observed DM energy 
density. The annihilation into two photons is smaller by a factor of about ten. The largest contribution to 
the annihilation cross section would come from either annihilation into Higgses or fermion-KK fcrmion pairs as 
discussed in Section [2] thus justifying the relic density computation discussed there. 
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(M$/A) 4 . For 1 TeV photons, HESS or VERITAS would be sensitive to A ~ 3.5 TeV. We 
summarize these observations in Table [U 

Ground-based Cherenkov detectors capable of operating at larger AQ can start probing 
theoretically interesting values of A even for halo profiles not as peaked as the Moore et. al. 
profile. From Fig. [3] we see that for AQ ~ 1, several halo profiles converge to J(AQ) x AQ pa 
10. An additional factor of 6 improvement in the flux sensitivity would then make scales 
A ~ 10 TeV accessible. Of course, larger AQ means also larger background, but hopefully the 
very characteristic line signal can be extracted if there are enough events (see Ref. |17j). 

We finally mention the possibility of observing photons from Higgs decays (assuming that 
the main channel for DM annihilation is into Higgses, as in Subsection 12.11 so that (cthhv) ~ 
0.8 pb). For instance, for a SM-like Higgs with mass around = 135 GeV, the branching 
fractions into 77 or Z7 are of order 10 -3 each. The photons from these channels present a flat 
spectrum between EffT, = |M$(1 — 0) and E^ x = |M$(1 + 0) for the 77 signal, or between 
E^ n = \M^{l-m 2 z /m 2 H ){l-0) and Ejg* = \M${l-m%/m 2 H ){l + l3) for the Z 7 signal. Here 
(3 = -i/l — mfj I m| is the velocity of the Higgs in the DM rest frame. 

We show in Fig. H] the total flux integrated from a threshold energy E^ = 50 GeV up 
to -Emax' as a function of M$. Here we optimistically assume J(Afl) x AQ = 10 2 as would 
be appropriate for the Moore profile though of course with other profiles the signal would be 
smaller. ACTs such as HESS or VERITAS would be sensitive to such a signal, but if the halo 
profile is less peaked or if the Higgs branching fraction into photons is smaller, this continuous 
signal becomes challenging. Nonetheless since this is a generic prediction of this type of model 
that doesnt rely on higher- dimension operators exploring the possibility of detecting such a 
signal is extremely worthwhile. 

3.2 Positrons 

Recently there has been intriguing evidence for an excess positron signal at energies up to about 
80 GeV j27j. Clearly it is of interest to determine whether such positrons can arise from dark 
matter annihilation. We do not anticipate that heavy dark matter particles will explain this 
excess, since a positron signal, if it exists, will be concentrated at higher energies. 

Nonetheless it is of interest to explore this positron signal to see how it compares to back- 
ground and to see whether in principle the signal could be detectable at high energies. We 
also briefly consider lighter particles (with less theoretical motivation in our context of strongly 
intereacting TeV scale physics) with similar interactions to KK dark matter particles and find 



16 



10 



10 




(cr HH v/c) = 0.8 pb 



10 



14 



200 



500 



1000 



2000 



5000 



m$[GeV] 



Figure 4: Integrated photon flux from — > HH and H — > 77 or — ► Z7, as a function 
of M$. We assume that the branching fractions for these decay modes are 10~ 3 and take 

J (An) xAn = io 2 . 

that even without a big boost factor annihilations of ~ 100 GeV dark matter matches the 
PAMELA data. 

With this in mind, we note that for the bulk scalar fields discussed in Subsection 12.21 other 
operators involving SM fields can be interesting from the point of view of DM signals, besides 
the operators leading to direct annihilation of the $ particles into photons discussed in the 
previous subsection. Specifically, the higher-dimension operators of the type Eq. (|Sj), coupling 
a pair of $s to an electron and its lowest KK mode, can lead to an interesting positron signal. 
The dominant annihilation channel involves the fermions closest to the IR brane (we discussed 
in Subsection 12.21 how the DM relic density can be determined by the annihilation into a SM 
fermion and one of its KK modes). 

In Subsection 12.21 we defined three scenarios that differ on how the fermions are localized in 
the extra dimension. Of these, the most favorable one to obtain a sizable positron signal from 
DM annihilation is scenario 3. But we will see that only the electrons need to be somewhat 
localized near the IR brane for the positron signal to be interesting, and this can occur in 
scenario 2 as wellfl Note also that the positron signal is sensitive to the local DM distribution 

9 A hybrid case (of scenarios 1 and 2) with one of the lepton chiralities localized somewhat near the IR 
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(and not very much to how peaked the DM halo is at the galactic center). It is common 
to parameterize the effects of DM inhomogeneities by an (energy-independent) "boost" factor 
B = (p 2 )/(p) 2 - Studies of such enhancements for gamma rays indicate that the boost factor 
might be as large as order 10 [28J. Therefore, the positron signal can receive an enhancement 
compared to the case of a smooth DM density distribution, though the likely size of this 
enhancement is not expected to be very large. 

With this understanding we proceed to estimate the signal from direct annihilation into a 
positron and a KK mode. The produced positron has a well-defined energy that depends on 
the DM and KK fermion masses, M$ and M eW : E^+ m ' = (AMI ~ M lm)/^ M $>- Besides these 
primary monoenergetic positrons, we also consider the secondary positrons arising from the 
annihilation of $ particles into an electron (neutrino) and a positron KK mode followed by the 
decay of the associated KK lepton into a positron and a Z (W) gauge bosonj^j When the KK 
lepton l^ 1 ' = or has a mass slightly below 2M$ (about the threshold for DM annihilation 
into lepton and KK lepton), it is produced nearly at rest and the resulting positron from its 
decay has a relatively well-defined energy. In detail, the energy of the KK lepton produced in 
DM annihilation is E lW = (4M| + M^ 1) )/4M $ , while its momentum is p = (4M| — Mf {1) ) /4M$ . 
For the two-body decays e^ 1 ' — > Ze or — > We, one finds the typical flat spectral distribution 

h(Eo) f E <*<E + (ig) 

[ otherwise 

where E is the positron energy. Neglecting the masses of the decay products, the maximum 
positron energy is E + = ^E t (i)(l + = M$, while the minimum positron energy is E_ = 
§^ C1) (1-Aci)) = M^/AM*. Here /%d =p/E m = (4M| - M^)/(4M| + Mf (1) ) is the velocity 
of the KK lepton (in the rest frame of $). Notice that the upper endpoint is determined by 



brane and the opposite chirality localized near the UV brane (to generate the small lepton masses by the 
exponential wavefunction suppression) also falls in this category. For instance if c; B f=s 0.4 while c ta ~ 0, we 
have T]t R /r]i R « [(1 — 2ct R )/(l — 2c/,,)] 1 / 2 ~ 2, where the 77 parameters were defined in Eq. Therefore, the 
annihilation into positrons can plausibly be suppressed compared to the dominant top channel by only a factor 
~ AN C ~ 10, if the unknown dimensionless coefficients A e and At in Eq. ([7]) are assumed to be comparable. 

If only the RH top and the RH leptons are localized near the IR brane and one neglects other annihilation 
channels, the thermal relic density computation implies (cr e (i) e v/c) = (cr„(i)„i?/c) = (a T (i) T v/c) ~ 0.06 pb and 
(o- 4 (D^/c) ~ 0.6 pb. 

10 The decays of the KK lepton into Higgs are suppressed by the electron Yukawa coupling. For gauge KK 
masses of order 3 TeV, the main decay channels of the KK lepton involve Z or W (through EWSB mixing of the 
Z/W with its KK modes, as opposed to mixing of the lepton and its KK modes). When the lepton is an SU(2) 
doublet we have T(eW -> Ze)/T{u^ -> We) » (r7i|/m^)(r 3 - s 2 w Q) 2 /c 2 w , leading to BR^ 1 ) -> We) « 75% 
and BR(e' 1 ' — > Ze) w 25%. Similarly, the SU(2) singlet KK positron decays dominantly into Ze + . 
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the DM mass only, and that in the limit M^i) — ► 2M$ one has ^(-^o) ~~ > — M$). Further 
decays of the Ws and Zs lead to additional positrons that have a softer spectrum and give a 
subdominant contribution due to the small branching fractions involved. We do not include 
positrons from W or Z decay in the following analysis. Note also that primary and secondary 
electrons with the exact same characteristics as the positrons above are also produced. 

The positron energy is distorted as it propagates through the interstellar medium be- 
fore detection. In general, for an initial spectral distribution fi(E ), normalized according 
to J °° dE fi(E ) = 1, the differential positron flux at the solar position is obtained from 



where po is the average DM mass density, B is the boost factor, (<7jt>) is the z-th channel 
thermally averaged (ultra non-relativistic) DM annihilation cross section times relative velocity, 
B l e+ is the corresponding branching fraction into positrons, and G(E , E) is a Green function 
that includes the details of the DM mass distribution in the galactic halo, takes into account 
the propagation of the positrons through the interstellar medium in the galaxy, and describes 
how their energy E is shifted under diffusion, various spatially and energy-dependent energy 
loss mechanisms, reacceleration, etc. The direct annihilation into positrons plus their lightest 
KK mode simply corresponds to f{E ) = 5(E — i^" m '), while secondary positrons arising from 
the decay of the KK lepton are described by Eq. ffT8l) . 

In Ref. [29], Moskalenko and Strong modeled the propagation of positrons through the 
interstellar medium for several galactic halo DM mass distributions. They provided a simple 
parameterization for the Green function that reproduces the more detailed simulation^"*] to 
within 10%: 



where G(Eq, E) is given in units of cm sr _1 GeV _1 , E is the local positron energy in GeV, and 
the coefficients a, b, c, w, x and y are functions of Eq (the initial positron energy) that are 
tabulated in Tables II and III of Ref. [29]. For defmiteness, we consider the "isothermal" model, 
which is characterized by a spherically symmetric DM mass distribution given by 



n The code used in the simulation aims at reproducing simultaneously observational data related to cosmic 
ray origins and propagation such as: direct measurements of nuclei, antiprotons, electrons and positrons, as 
well as indirect measurements via 7 rays and synchrotron radiation. 




(19) 



P(r) = Po 




2 ' 



(21) 
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where r c is the core radius and R Q = 8.5 kpc is the solar distance to the galactic center (the 
parameters r c and po are obtained by fitting to the rotation curve, and for the isothermal model 
r c = 2.8 kpc and po — 0.43 GeV cm" 3 ). We use a galactic halo size of Zh = 10 kpc, which is on 
the upper limit of the 4 — 10 kpc range favored by the analysis in [29] . 
The local positron flux then takes the form 

where the units on the r.h.s. are GeV cm -2 s~ 1 sr~ 1 , and the dimensionless Fi(E) is defined by 

Fi(E) = 10~ 25 E 2 J dE f t ,(E )G(E , E) . (23) 

Experimental observations are commonly presented in the form of the positron fraction, 
e + / (e~ + e + ), where the electron and positron fluxes include both background and signal. This 
quantity has the advantage that systematic uncertainties cancel out. For the electron and 
positron background spectral distributions we use the simple parameterizations given in [30] : 

/ d<S> e - \ 0.16 £- L1 

v^^; prim . bkg 7_ i + n^ 9 +3.2^.i5 > 

/ d$ e - \ 0.70 E - 7 

\dQdEj sec bkg ~ 1 + 1 10 E 1 - 5 + 600 E 2 - 9 + 580 E*- 2 ' 

(d$ e+ \ 4.5 E - 7 

\dQdEJ sec bkg ~ 1 + 650 E 2 - 3 + 1500 E 4 - 2 ' 

where E is in GeV and the units of the l.h.s are GeV -1 cm" 2 s" 1 sr -1 . 

In the following we envision a scenario in which the dominant DM annihilation channel is 
into pairs of electron/positron plus a KK mode, as could be expected in scenario 3 defined in 
Subsection 12.21 In this case, the WMAP relic abundance requires {v e w e v/c) ~ 1 pb. Somewhat 
more generally, the results are valid for Bx (a e (i) e v/c) = 1 pb, where B is the boost factor. Since 
B is expected to be order a few, our results can illustrate situations where the electron/positron 
channel is one among a few dominant annihilation channels (e.g. if the other lepton channels 
are equally important). 

In the left panel of Fig. [5] we show the expected positron fraction signal, assuming DM 
masses M$ = 500 GeV and M$ = 2 TeV. We also include the secondary positrons arising 
from the decay of the KK lepton, assuming M e (i) = M$. Interestingly, there is a rather clear 
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Figure 5: Left panel: positron fraction including the primary positrons/electrons from the 
annihilation $$ — > e ± e^ and the secondary positrons/electrons from the annihilation $$ — > 
/W/( ', followed by a two-body decay — > e ± X. We show the spectra for two DM masses, 
M$ = 500 GeV and M$ = 2 TeV, assuming that M e (i) = M$ and a boost factor such that 
B x (a e (i) e v/c) = 1 pb, with p = 0.3 GeV/cm 3 . The solid lines represent the individual 
contributions from primary and secondary production. Right panel: positron fraction for M$ = 
1 TeV and M lW = 0.85 x (2M$) = 1.7 TeV showing clear peaks at E^ m - = (4M| -M e 2 (1) )/4M$ 
and near M$ (the endpoint is exactly at M$). 

peak above background that could be observable in the sub TeV range. A moderate boost 
factor would make such a feature even more prominent. In the right panel of Fig. [5] we show 
another example where M e (i) = 0.85 x 2M$ is closer to threshold. This case exhibits more 
clearly the two peaks discussed above, one at E w ™' from the primary positrons and another 
near M<j, for the secondary positrons from the KK lepton decay. As remarked above, the upper 
endpoint gives a direct measurement of M$. The first peak then gives information about the 
KK lepton mass. With a handle on the DM and KK fermion masses it would be possible to get 
information about the effective cutoff scale (A e /AL) _1 A [see coefficient a in Eq. (TTTj) ] . modulo 
the uncertainty associated with the local DM energy density (since these high-energy positrons 
come from distances of at most a few kpc, the dependence on the DM halo model is expected 
to be milder). 

There seems to be evidence in several experiments for an excess in the positron flux in the 
tens of GeV energy range (HEAT {31], AMS-01 [32]), with the PAMELA satellite experiment 
supporting this excess up to energies of about 80 GeV [27J . As has been emphasized recently [33, 
[341, 135] the observed fluxes are larger than what would be expected from thermal WIMPs when 
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Figure 6: Positron fraction due to positrons from a 120 GeV DM particle annihilating into ee^\ 
where eW is a heavy electron of mass M e (i) = 200 GeV. The lower energy peak arises from 
direct annihilation into a positron (plus heavy lepton), while the higher-energy peak arises from 
the direct decay of the heavy lepton into a positron (plus a gauge boson). Further positrons 
from the W's or Z's produced in the heavy lepton decay are not included and would add a 
softer contribution. We assume B x (<r e m e v/c) = 3 pb. 



these lead to positrons mainly through the decays of their annihilation products (e.g. W's 




Nonetheless, it is interesting to note that a DM candidate with a mass of about 100 GeV 
annihilating primarily into electrons/positrons can explain the observed positron excess with a 
boost factor of order unity. Although such low masses are not expected in our scenario, we show 
in Fig. [6] the positron signal from the annihilation of a 120 GeV DM candidate into a positron 
and a "heavy vector-like electron" of mass M e (i) = 200 GeV. We show the HEAT data and 
the recently released PAMELA data, which shows a clear increase with energy of the positron 
fraction up to energies of at least 80 GeV [27]. However, the ATIC-2 balloon experiment [38J 
also indicates an excess in the total electron plus positron flux extending up to energies of about 
a TeV, which would not be explained by the self-annihilation into electrons/positrons of such 
a light DM candidate. Antiprotons produced in decays of the heavy electron (via W gauge 
bosons) may also conflict with the non-observation of an antiproton excess in the PAMELA 
p/p data [39] . 

Of course the above observations hold only for a light dark matter particle of order 100 GeV. 

12 However, the presence of a relatively long range force among the DM particles can lead to an enhancement 
of several orders of magnitude in the annihilation cross section at very low velocities, which can account for these 
observations in certain dark matter models [341 136 [ [37j. Our DM candidate, having only non-renormalizable 
interactions, does not fall into this category. 
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Figure 7: Positron fraction when —>■ WH (including longitudinal gauge bosons) dominates 
the DM annihilation cross section. We include the secondary positrons/electrons from the 
annihilation $$ — ► W£\V£ and the tertiary positrons/electrons from the annihilation — ► 
hh — > 4W, assuming nih = 170 GeV. Only two-body W decays into positrons/electrons are 
included. We take M$ = 1 TeV and B x (a w+w -v/c) = 10 pb, with p = 0.3 GeV/cm 3 . The 
solid lines represent the individual contributions from secondary and tertiary production, with 
the latter having a softer spectrum. 

For the range of masses expected in our scenario (around 1 TeV), the proton flux is sufficiently 
suppressed, but the positron excess below about 80 GeV would not be explained. However, the 
signal can exceed the background at high energy which would be interesting if experiments can 
attain the required sensitivity. 

We close this subsection by coming back to the case of a scalar DM candidate that decays 
dominantly into Higgses, as discussed in Subsection l2.11 and comment on the associated positron 
signal. The annihilation cross section Eq. ([3]) corresponds to the processes — >• WlWl, 
$$ —> Z L Z L and — > hh in the large M$ limit. In this limit we have a(§§ —> W L W L ) pa 
2cr($$ — > Z^Zl) m 2cr($$ — > hh). Further decays of the Ws, Zs and Higgses can result in 
energetic positrons. We consider here the case m# = 170 GeV with a SM branching fraction 
BR(if — > Vy + Vy~) ~ 1, and compute the spectrum of secondary positrons from — > 2W/Z — > 
e + X and of tertiary positrons from $<3> — * hh —> 4W — > e + X. We do not include positrons 
from processes further down the decay chain. 

The cross sections times branching fractions (including the positron multiplicities) to be used 
in Eq. ([22]) are (a w+w -v) x BR(W -> eu), 2 x (a zz v) x BR(Z -> e+e") w (a w+w -v) x BR(Z -> 
e+e"), and 2 x (cx^) x BR(H -> W+W - ) x BR(iy -> ez/) w (avi/+iy-w) x BR(W ei/), 
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respectively, where BR(W — > eu) ~ 0.11 and BR(Z — > e + e~) ~ 0.036. In Fig. [7J we show 
the positron fraction for M$ = 1 TeV and B x {a w + w -v / c) = 10 pb, which corresponds to a 
boost factor B w 20. We conclude that such a signal would be visible above background only 
for rather large boost factors. In particular, such a scenario can also not explain the observed 
positron excess reported at lower energies. 

4 Conclusions 

We considered a simple scenario for scalar DM of mass around 1 TeV in the context of non- 
renormalizable theories with a cutoff near the TeV scale. Such a possibility arises naturally in 
extra-dimensional models that address the hierarchy problem, such as the Randall-Sundrum 
scenario but could also arise from other strongly interacting TeV-scale theories. The thermal 
relic density can be determined either by renormalizable or non-renormalizable interactions. 
Such dark matter particles are clearly more challenging to detect but can conceivably yield ob- 
servable gamma ray signals at current detectors and might ultimately yield observable positrons. 

A monochromatic gamma ray line signal arises from the direct annihilation via nonrenor- 
malizable operators of the DM particle $ into photons. For a cutoff scale of up to about 10 TeV, 
such a signal can be larger than the signal from a typical one-loop induced direct coupling to 
photons. We also point out that the monochromatic signal associated with non-renormalizable 
operators is likely observable in currently operating ground-based experiments, and could be 
used to probe the cutoff scale up to several TeVs. 

We also consider secondary photons from Higgs decay in the annihilation — > HH, and 
find that this continuous signal can be observable if the DM halo is relatively peaked at the 
galactic center. Secondary or tertiary positrons can also be produced in the decays of Higgses or 
longitudinally polarized gauge bosons, but the positron flux is likely too small to be observable 
above background. 

It is also possible to have more exotic scenarios where the annihilation cross section is 
dominated by the non-renormalizable interactions, as opposed to the dimension-4 coupling of 
$ to the SM Higgs field. In the extra-dimensional context one could have annihilations into a 
SM fermion and the associated KK fermion dominating the total annihilation cross section. If 
the leptonic channels are dominant, it is possible to have an observable positron signal in the 100 
GeV to 1 TeV range with a boost factor of order one. Such a signal would typically present two 
peaks, due to the heavy lepton involved. However, the expected flux is too small to account 
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for the positron excess reported by HESS/AMS-01 and the PAMELA satellite experiment. 
Nevertheless, we find it promising that indirect searches in the sub-TeV range can be sensitive 
to cutoff scale physics. 
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A Bulk Scalars in an RS Background 

We consider a bulk real scalar, $, propagating in the background j5] 

ds 2 = e- 2ky Vflu dx^dx u - dy 2 , (24) 

where x M (/i = 0, 1, 2, 3) are the 4D coordinates, and < y < L parametrizes the fifth dimension. 
We assume that the scalar obeys (— , +) boundary conditions, and consider the action 

S = J d 4 x J dy^/\i\^ {d M $d M <& - M 2 $ 2 - S(y - L)m$ 2 } , (25) 

where M and m are bulk and IR localized mass parameters, respectively. We do not write 
a localized mass on the UV brane, since the scalar is assumed to vanish at y = 0. We will 
parametrize these mass parameters in units of the curvature scale as 

15" 
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rn 



c 2 s + c s - i 

3 r 



k 2 , (26) 
k . (27) 



For 5 = 0, the mass spectrum that follows coincides precisely with that of a fermion obeying 
(— ,+) b.c, where Cf = c s parametrizes the fermion bulk mass [40J. Our sign conventions are 
such that for c s < 1/2 the lightest KK mode is exponentially localized near the IR brane, 
and when c s < —1/2 its mass is exponentially smaller than the warped down curvature scale 
k = ke~ hL . We will see that the lightest eigenvalue can remain small for a wide range of values 
of the parameter 5 defined in Eq. (1271) . and therefore the lightest scalar KK mode can be easily 
lighter than the SM gauge and fermion KK resonances. 
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The KK decomposition for $ reads 

ky ~ 

^ i y) = -=Y i ^)fM , (28) 

where we pulled out an explicit factor e ky for convenience, and the KK wavefunctions obey 

d% - 2kd y f n - (3k 2 + M 2 )f n = -e 2ky m 2 J n , (29) 

and satisfy the b.c: 

f n (0)\ y=0 = , d y f n \ y =L = ~(k + m)f n (L) . (30) 
The solutions can be written in terms of Bessel functions as 

fn(y) = A n e*y [j |Cs+|| (^) +6^, (!£e*)] , (31) 
where A n is a normalization constant, determined from 

1 r L 

and 



J dyfn(y)fm(y) = S nm , (32) 



J \c + 1 l( I T t ) 

&=- |a+ *' * . (33) 

r |c s + i|l fc ) 



The eigenvalues can be written ke hL , where the x n solve 

J\ Cs+ i\(x n e- kL ) x n J ]Cs+ i hl (x n ) + (c s + | - \c s + || +5)J| Cfl+ i|( 



,1 ■■■■>■, 



^ Cs+ i|(s n e feL ) x n ^ Ca+ i hl (x n ) + (c a + | - \c s + || + 5)^ Cs+ i|(x n 



(34) 



For E>1, the lowest solutions are approximately given by the vanishing of the numerator 
in the r.h.s of Eq. (j34H . We show the smallest eigenvalue in Fig. [8] as a function of 5, defined in 
Eq. ([27)1 . for several values of c s [which parametrizes the bulk mass M 2 as in Eq. (|26|) ]. This 
eigenvalue can be well approximated by 



■<>M V V ( 35 ) 

2 vf + c * v c s > -1/2 . 
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Figure 8: Lightest mass, mi = xike~ kL , for a bulk scalar obeying (— , +) b.c., as a function of 
bulk and IR localized masses, parametrized by c s and 5, as in Eqs. ( 1261) and ( 1271) . Recall that 
for a KK gauge boson obeying (+, +) b.c, x\ ~ 2.45. 

We see that for 5 < Min{0, —(1 + 2c s )}, x\ becomes negative and the corresponding mode is a 
tachyon. We will assume that we are in a region where no such instability arises. One should 
also keep in mind that for 5 = and c s ~ —1/2, Eq. (135j) receives additional corrections not 
shown there. In this case, the smallest eigenvalue remains non-zero, becoming exponentially 
small for c s < —1/2. 

For Fig. [2]in the main text, we chose c s = —0.2 and adjusted S so as to reproduce the desired 
mass M$. This determines the corresponding wavefunction and allows the computation of the 
relevant overlap integrals that determine the $ couplings. Note, however, that the dependence 
on the choice c s = —0.2 is very mild. 
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